Nanorheology of viscoelastic shells: Applications to viral capsids 
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We study the microrheology of nanoparticle shells [Dinsmore et al. Science 298, 1006 (2002)] 
and viral capsids [Ivanovska et al. PNAS 101, 7600 (2004)] by computing the mechanical response 
function and thermal fluctuation spectrum of a viscoelastic spherical shell that is permeable to 
the surrounding solvent. We determine analytically the damped dynamics of the shear, bend, 
and compression modes of the shell coupled to the solvent both inside and outside the sphere in 
the zero Reynolds number limit. We identify fundamental length and time scales in the system, 
and compute the thermal correlation function of displacements of antipodal points on the sphere 
and the mechanical response to pinching forces applied at these points. We describe how such a 
frequency-dependent antipodal correlation and/or response function, which should be measurable 
in new AFM-based microrheology experiments, can probe the viscoelasticity of these synthetic and 
biological shells constructed of nanoparticles. 

PACS numbers: 87.15.La, 61.46.-w, 82.70.Uv 
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I. INTRODUCTION 

Understanding the dynamics of soft, nanoporous ma- 
terials will be an important area of research both in 
biophysics and in the materials science communities. 
This emerging field of study is being driven by new 
AFM-based mechanical or rheological measurements of 
supramolecular biological materials [l], 0] and even entire 
cells p| H| . Perhaps the prototypical examples of such a 
material are found in the plethora of viral capsids. Such 
shells are constructed from a small set of proteins or their 
oligomers (capsomeres) arranged in an ordered structure 
forming roughly a spherical shell. The packing of these 
capsomeres generically leaves at least nanometer scale 
pores through which solvent may flow. In addition, there 
is a wide variety of similarly nanoporous synthetic struc- 
tures such as linked networks of nanoparticles Q and 
colloids @, 0] assembled into two-dimensional shells. In 
both the biological and synthetic cases, the material mak- 
ing up the thin shell or membrane is constructed from 
identical nanoscale tiles that incorporate holes of com- 
parable size through which the solvent may flow. Other 
examples of porous membranes include lipid bilayers con- 
taining pore-forming transmembrane proteins [1, Q and 
even the large (~ lOOnm) cytoplasmic ribonucleoprotein 
vaults [l(| , which make up a rather ubiquitous but enig- 
matic component of eukaryotic cells. 

We focus on elucidating the mechanical response of 
such a permeable shell to a set of sinusoidally time- 
varying "pinching" forces applied at antipodal points on 
the sphere - see Fig. [IJa) . Such forces represent the sim- 
plest characterization of an AFM-based nanoindentation 
experiment performed at finite frequency. By applying 
known forces at a fixed frequency and observing the in- 
phase and out-of-phase response of the diameter D of the 
spherical shell one should be able to extract the viscoelas- 
tic properties of the shell material. Such an experiment is 



the finite-frequency generalization of the work of Michel 
et al. [Tl| and can be considered to be an active microrhe- 
ological measurement. Alternatively, one may imagine 
that the AFM tip can be used in a passive manner to 
monitor the thermal fluctuations of this diameter. Since 
a number of these physically relevant shells are rather 
incompliant, their thermal fluctuations are likely to be 
unresolvable. In such cases passive microrheology (i.e. 
the monitoring of thermally fluctuating variables) must 
be supplanted by active techniques. This distinction is 
of little importance for our analysis provided that the 
mechanical perturbation in the active technique is small 
enough that the system remains in its linear response 
regime. In the passive technique one has direct access 
to the power spectrum of the fluctuating quantity deter- 
mined (via the fluctuation dissipation theorem) by the 
imaginary part of the response function. The real part is 
then obtained by a Kramers-Kronig relation [l§] ■ In the 
active measurement the real and imaginary parts of the 
response function are found directly from the in-phase 
and out-of-phase response of the shell respectively. In 
this article we will compute the full frequency-dependent, 
complex response function as well as the expected ther- 
mal power spectrum of (\D(u})\ 2 ) of each shell studied. 

To understand the mechanics/dynamics of such 
nanoscale objects in solution, one must consider the cou- 
pled dynamics of soft membranes/shells and the sur- 
rounding solvent. The case of viscous shells that are 
impermeable to the surrounding solvent has been well- 
studied in the context of microemulsions and vesicles [l3| . 
The focus of our current study, however, is the dynamics 
of membranes or shells that are permeable to the sur- 
rounding solvent on time-scales relevant to the experi- 
ment. The effect of solvent permeation, which we quan- 
tify below, is to create a new dissipative stress on the 
shell related to solvent permeation. These dissipative 
(out-of-phase) stresses associated with fluid permeation 
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FIG. 1: (color online) a) A schematic illustration of an AFM- 
based rheological experiment on a nanoscale shell. The (red) 
arrows depict the applied stresses ni = —772 that deform 
the (blue) shell of radius R from a sphere on the left into 
an oblate object shown in the right. The total change in 
sphere's diameter along the line of force application is SD. 
This deformation can be resolved into a linear superposition 
of two sets of normal modes. In b) we show the Lennon- 
Brochard t — 2 mode while in c) we show the compression 
dominated longitudinal sound £ — 2 mode. In each figure the 
(black) arrows represent the velocity field of the sphere. 



must also be taken into account for the quantitatively 
correct interpretation of such measurements as described 
above. In particular in the case of a membrane or shell 
that has an inherent viscoelastic response, one may ask 
whether it remains possible to distinguish the dissipative 
stresses due to from the material from those associated 
with fluid permeation. 

In this article we quantify and explore the effects of 
membrane permeability on the dynamics of both flat and 
spherical membranes. We parameterize the continuum 
mechanics of the shells by a permeation coefficient, a 
bending modulus, an area compression modulus, and a 
two-dimensional, in-plane shear modulus. These latter 
two moduli we will allow to be viscoelastic, i.e. com- 
plex and frequency dependent. We discuss the interpre- 
tation of this generalization below. We show that solvent 
permeability leads to a new dissipative stress in the sys- 
tem that changes the relaxation rates of each deformation 
mode of the system. We show also that, on a spherical 
shell, these normal modes are superpositions of compres- 
sion and bending deformations. We then consider four 
test cases to explore the interrelated roles of shell poros- 
ity and viscoelasticity in determining their mechanical 
response. 

First, we examine the expected response of a purely 
elastic, porous shell to isolate the dissipative effects of 
the interaction of the porous shell with the surrounding 
solvent. We then consider the case of a highly incom- 
pressible, but purely viscous membrane of varying poros- 
ity. Such systems are reminiscent of giant unilamellar 
vesicles (GUVs) [0], [HI, EH containing pore- forming pro- 
teins [13, HI, Hal - 



Finally, we consider two simple models of a viscoelastic 
shell inspired by viral capsids. In the first case we assume 
that the shell has a purely elastic response to compres- 
sion but a viscoelastic response to in-plane shear stresses. 
Such a viscoelastic response to in-plane shear may result 
from plastic rearrangements of capsomere proteins under 
applied stresses. Thus, the proteinaceous shell of a virus 
may deform like an amorphous solid. In the second case 
we assume an elastic response to in-plane shear stress, 
but the dissipative relaxation of in-plane compressional 
stress. Such stress relaxation should occur in cases where 
the capsomeres have internal degrees of freedom associ- 
ated with allosteric conformational transitions [2fl l2~il. 1221] 
that change their cross-sectional area or packing density. 
To model this effect we assume a compressional stress re- 
laxation time on a scale typical of known allosteric tran- 
sitions in proteins. In both cases we employ the simplest 
model of viscoelasticity, the Maxwell model (23[, which 
has a single stress-relaxation time scale. 

The central question we address is whether it is pos- 
sible to extract the dissipative mechanics of the capsid 
proteins from the proposed AFM-based microrheological 
experiment. It is possible that the dissipative stresses as- 
sociated with internal rearrangements of the capsomeres 
or plastic deformations of the shell will be lost against the 
background of dissipative forces associated with solvent 
flow through the porous membrane. In both cases, owing 
to a separation of time scales, the nanoindentation exper- 
iment should be able to resolve the expected viscoelastic 
response of the shell. We show, however, that in the for- 
mer case of a viscoelastic shear response, the mechanical 
signature is quite subtle in the real part of the response 
function due to the dominance of the bending and com- 
pression moduli in the response function. However, the 
imaginary or out-of-phase part of the response function 
is capable of recording this effect. In the latter case the 
effect of compressional stress relaxation on time scales 
typical of protein conformational change should give a 
dramatic mechanical signature in the imaginary part of 
the response function so that finite-frequency nanoinden- 
tation studies should be an excellent probe of these com- 
plex mechanics. 

The remainder of the paper is organized as follows. 
In section Hi Al we calculate the effect of porosity on the 
undulation dynamics of a flat membrane in a viscous sol- 
vent. We also discuss expected values of the permeation 
coefficient for a variety of physically important systems. 
The calculation of these coefficients is discussed in ap- 
pendix [X] In section III Bl we consider the deformations 
of a spherical membrane where the curvature of the un- 
deformed structure couples the previously studied bend- 
ing mode to in-plane compression. To incorporate the 
mechanical effects of solvent hydrodynamics in the case 
of spherical shells, we recount and apply the pioneer- 
ing work of Lamb [2J] and Brenner [25| on low-Reynolds 
number flows inside and outside a sphere in section [II CI 
Using the results of the previous calculations we plot and 
discuss the decay rates of the linearly independent modes 
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of the shell in section III Dl In section IIIII we use the dy- 
namics of these linearly-independent modes to construct 
the complex, frequency-dependent response function of 
the shell to antipodal pinching forces. In sections IlII Al 
IIIIB1 and IIII Cl we plot the response functions and pre- 
dicted thermal power spectra of elastic, viscous, and the 
two classes of viscoelastic shells introduced above. Fi- 
nally, in section HVl we compare these results and discuss 
the implications of our work for future nanoindentation- 
bascd rhcology measurements. 



II. CALCULATION 
A. The plane permeable membrane 

The overarching feature of these systems is that fluid 
flow through the shell or membrane (permeation) can oc- 
cur on timescales comparable to the relaxation of surface 
deformations. Thus, the velocity of the fluid at the mem- 
brane surface in the direction along the local normal of 
that surface will not be equal to the normal velocity of the 
membrane. In the limit of high membrane porosity, the 
fluid will essentially pass though the membrane undis- 
turbed. Consequently the hydrodynamic interaction of 
the membrane with itself, which is known to qualitatively 
alter the spectrum of the decay rates of membrane un- 
dulations [26| is suppressed and replaced by dynamics 
consistent with local drag. As the porosity or permeabil- 
ity of the membrane is reduced, the system continuously 
recovers these well-known hydrodynamic self-interactions 
that control the structure of the spectrum of undulatory 
decay rates. To better understand this feature we ad- 
dress the dynamics of a plane porous elastic membrane 
in a viscous Newtonian solvent. We assume that the dy- 
namics are all occurring at low Reynolds number. For 
a flat membrane or interface the normal modes of the 
coupled interface and solvent system involve bending in- 
plane shear and compression [27| . Since we will concen- 
trate on indentation-based measurements, we consider 
only the undulation mode of the flat membrane. 

A zero-tension membrane is characterized by a single 
bending modulus k so that free energy of the membrane 
may be written as 



' dxdy, 



(1) 



in terms of normal displacement w. Here is the in- 
plane gradient operator acting on membrane coordinates 
x,y. In the above we implicitly assumed small deforma- 
tions of the surface from its flat equilibrium shape. The 
restoring stresses acting on the deformed membrane are 
given by 
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(2) 



The equation of motion for the fluid of viscosity rj in 
the limit of vanishing Reynolds number is the linearized 



Navier-Stokes equation for an incompressible fluid: 
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(3) 
(4) 



where p is the pressure. The interaction between the fluid 
and the membrane is governed by three boundary con- 
ditions enforced at the membrane lying in the xy plane. 
First, we require Darcy's law to be obeyed so that the 
flux of the fluid through a unit area of the membrane is 
proportional to the pressure difference between the two 
sides [28]. Thus 



dw 
~dt 



P 



P 



(5) 



where p + '~ are the pressures above (z > 0) and below 
(z < 0) the membrane respectively. We discuss the es- 
timation of the Darcy coefficient a in terms of the mi- 
croscale structure of the membrane in Appendix [XJ Ex- 
pected values of the permeation coefficient for a variety 
of systems are given in Table 1. This first boundary con- 
dition relates the normal velocities of the membrane and 
fluid at the membrane surface by equating the normal 
stress difference across the membrane to the fluid flow 
through it. In addition we must insist upon the conti- 
nuity of the normal component of the fluid velocity so 
that 



v+(Q)=v-(0), 



(6) 



where vf(0) = lim 2 ^ ± v z (z). For the tangential veloci- 
ties, we enforce a no-slip boundary condition on the fluid 
at the upper and lower boundaries of the membrane: 



«+(o) = «-(o) = o. 



(7) 



Finally, we require the fluid velocity field and pressure to 
vanish infinitely far from the membrane so that 



v(±oo) = 0, p(±oo) = 0. 



(8) 



Ignoring the inertia of the membrane, we require 
force balance between the elastic restoring stress / from 
Eq. @, the viscous stresses cr+ and a~ z exerted by the 
fluid on both sides of the membrane and an externally 
imposed stress, ir e so that 



(Vzz - viz) 



(9) 



To compute the decay rate of membrane undulations 
having wave vector q, we consider a sinusoidally applied 
normal stress 7r e 



7r„ = n e tqx ~ lult , 



(10) 



and look for solutions of the membrane deformation tak- 
ing form of 



w 



w e 



iqx—iujt 



(ii) 
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TABLE I: Dimensionless permeation coefficients (Eq. 
for nanoscale capsules 



6.1) 



System 


Permeation coefficient 


Ref 


Rotavirus 


50 


[29] 


S-layer shell 


10 4 


[30] 


Colloidosome 


10 4 - 10 5 


m 


GUV 


10 - oo 


[8,9] 



Details of the calculation are presented in Appendix [B] 
In the limit that viscous stresses dominate inertial ones 
in the fluid, i.e. top <C rjq 2 , we find that the amplitude 
of the driven undulations is given by 



7T 


a 

2qr) ' 


f 2 




2 i oj a — 


Kq 4 


a 
2qr, 


1-2 



(12) 



w 



The poles of the response function shown in Eq. (|T2|) , de- 
termine the relaxation rate of the freely decaying surface. 
We find that the decay rate takes the form 



Si 



nq 
a 



(13) 



In the limit of an impermeable membrane, i.e. a — ^ oo, 
we recover the Lennon Brochard result [261 ] - As a be- 
comes smaller reflecting the increasing permeability of 
the membrane, the decay rate that scales as q 3 crosses 
over to q A scaling. This small-a dynamics is consistent 
with a "free-draining" assumption that the fluid merely 
exerts a local drag on the membrane, but does not con- 
tribute to a long-ranged hydrodynamic interaction be- 
tween distant parts of the surface. For any finite perme- 
ability, one may note from Eq. (| 13[) that there is a cross- 
over length I ~ rj/a above which the decay of surface 
undulations is controlled by long-range hydrodynamic in- 
teractions, while undulations of a wavelength less than I 
experience the free-draining dynamics. 

Given this intuition gained from the more simple ques- 
tion of small undulations of a planar, permeable mem- 
brane, we now turn to the problem of a spherical shell 
that is a more appropriate model of a viral capsid or its 
synthetic mimics. 



B. Mechanics of a porous spherical shell. 

Extending our analysis of the dynamics of a plane 
porous membrane to the case of spherical ones involves 
two separate complications. The first is due to the spheri- 
cal geometry of the undeformed membrane. In the case of 
a flat interface the in-plane shear and compression modes 
of the material decouple from the out-of-plane bending 
modes to linear order in the bending deformation [3l| . 
This is not the case for the sphere. For small deforma- 
tions of surfaces having finite curvature, out-of-plane de- 



formation is coupled at first order to in-plane dilatation 
and shear [32| . 

We write the elastic bending free energy JF fc of the 
spherical shell using the well-known Helfrich form [l3l . 
MM 



d 2 s[K a a - c ) 



(14) 



where d 2 s is an element of area on the surface, k is once 
again the bending rigidity, K% is the trace of the curva- 
ture tensor, and cq is the spontaneous curvature. The 
Greek indices run over the (angular) coordinates of the 
undeformed sphere. The curvature tensor on the spheri- 
cal membrane of radius R may be written in terms of the 
out-of-plane deformation, w, as 



K$ = D a Dpw + 5%w/R 2 



6? 3 /R. 



(15) 



In order to address out-of-plane deformation on the 
sphere, we must include the in-plane elastic stresses as- 
sociated with the dilatation and shear of the membrane. 
To do this we introduce a two-dimensional (2D) covariant 
strain tensor given by 

E% = ^(D a t p + D fj t a )+6%w/R, (16) 

where t is a displacement vector in tangent space of the 
sphere. Using this decomposition of the deformation 
into in-plane and out-of-plane motion, the most general 
displacement of material elements of sphere can be ex- 
pressed in the form 



w r + t. 



(17) 



The in-plane elastic energy of the sphere can be written 
as a surface integral of a scalar elastic energy density 



d 2 s 



X 



(18) 



where [i and A are the two 2D Lame coefficients required 
to describe the elasticity of isotropic materials. It is not 
clear a priori that the isotropic elasticity is sufficient to 
accurately describe the mechanics of a viral capsid or 
other ordered structures, but this choice introduces the 
minimal set of unknown parameters and greatly simpli- 
fies the analysis. In discussions of the equilibrium spec- 
trum of shape deformations of these membranes we will 
consider the viscoelastic generalization of Eq. (fT5)) where 
the Lame coefficients are complex functions of frequency 
so that the strain energy involves an integral over the 
history of prior deformations. 

Combining Eqs. (fT6| and (|18jl it is clear that the strain 
energy contains a term bilinear in out-of-plane and in- 
plane deformation. This shows that the curvature of the 
surface does indeed induce a linear order coupling of the 
in-plane deformation to forces along the local membrane 
normal. The total elastic free-energy J- c \ of membrane is 
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then expressed as the sum of Eqs. (fT4f and (fTSf . i.e. it 
is given by T c \ = T h + T e . 

We now look for the normal modes of the deformation 
of the sphere. Based on intuition from flat space, we 
write the in-plane displacement field t as the sum of an 
irrotational and a solcnoidal part 



tp = Dp* + la pD a X , 



(19) 



where ^ and x are two scalar fields defined on the surface 
of the sphere and 7 Q ^ is the alternating tensor. From 
Eq. (fT9)) one notes that \& determines the dilatational 
deformation of the system, while x describes the density- 
preserving shear modes of the membrane. 

Using this decomposition, the elastic free energy of the 
spherical shell is given by T e \ = Ti[w,ty] + ^[x]) where 

T 1 [w,^]= J d 2 s^w{A ± + 2/R 2 ) 2 w 

+ ^2 ( 2 / R -c^w + ^w 2 + —W A^ 

+ ~(K + » l )(A ± *) 2 + j^*(A x *)\ (20) 



and 

= J d 2 s {^A ±x [Axx + ?x/R 2 ] } ; (21) 

Aj_ is the two-dimensional in-plane Laplacian and K = 
[i + A 32]. Since there is no coupling between {w,^} 
and x, purely radial deformations of the sphere will ex- 
cite only the compression modes related to <!'. There- 
fore, in the interest of studying microrheological ap- 
proaches to the measurement of the membrane mechan- 
ics via nanoindentation studies, we may neglect ^[x]- 
Of course, to address the question of the dynamics of 
the porous, spherical membrane we must consider the 
(visco-) elastic object described above coupled to the sol- 
vent flows around and through it. In order to do so, we 
examine the fluid motions inside and outside the sphere 
in the zero Reynolds number or creeping flow regime. 



C. Hydrodynamics of the spherical shell 

It has long been recognized that studying this hydro- 
dynamics of a fluid confined to either the interior or ex- 
terior of a sphere is facilitated by recasting the Stokes 
equation in spherical coordinates and expanding the solu- 
tions in a basis of the appropriate solid spherical harmon- 
ics 0, [H, HIJ ■ These results were used 35] to explore 
the dynamics of an impermeable elastic shell immersed 
in a viscous fluid. In the interests of presenting a self- 
contained analysis, we recapitulate some of the earlier 
work on the hydrodynamics of the problem. We then dis- 
cuss and expand upon the work of Schneider et al. [35l[36j 
to introduce permeability, and to correct the mechanical 



coupling between bending and in-plane dilatation that is 
imposed by the curvature of the membrane as demon- 
strated in Eq. ([TH]) . 

To set the boundary conditions in a manner similar to 
that discussed for the case of an undulating permeable 
plane, we enforce the continuity of the radial components 
of the fluid velocity normal to the surface of the mem- 
brane so that 

<\ r=R =vr\ r=R , (22) 

where here and in the following we label the field asso- 
ciated with the fluid in the interior by "in" and on the 
exterior of the spherical shell by "out." We also require 
a no-slip boundary condition on the components of the 
velocity field in the tangent plane of the membrane so 
that 

v l l =vl ut = Vj_, (23) 
where the velocity of the membrane V is given by 

de dw „ _ dib . . 

Y = Tt = Tt r + v ^ (24) 

using Eq. (IT)) . Here Vj_ is the gradient operator in the 
tangent plane of the sphere and r is the local outward 
unit normal. 

To couple the flow of the fluid through the membrane 
to the deformation of it, we insist on force balance at the 
membrane so that 



— r fluid, 



(25) 



where F is the viscoelastic restoring force acting on the 
membrane due to its deformation history, and Ffl ui d is 
the negative of the difference between the viscous stresses 
m,out on ^ e mein |- )rane (j ue ^0 fluid flow inside and out- 

side of it respectively, 



(F fl 



uid)j = r k [a ko 



(26) 



Finally, the local difference in the normal velocity of the 
fluid and the membrane is governed by Darcy's law 



dw 
~dt 



(27) 



where a is the permeation parameter introduced in 
Eq. ©. 

The rotational symmetry of the problem ensures that 
the coupled membrane deformations and fluid motion 
corresponding to each spherical harmonic decouple from 
all other such terms. We can expand the fluid velocity 
field inside and outside the spherical membrane in terms 
of two scalar potentials, $, x f an d pressure p 0, [H, [HI]. 
Due to incompressibility and the Stokes equation, one 
can show that these functions and the pressure field p 
each satisfy a Laplace equation, 



V 2 p = 0, V 2 $ = 0, 



vY 



o. 



(28) 
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It can be shown that the field x f is the source of fluid vor- 
ticity. To analyze the effects of indentation experiments 
on the sphere, we will study membrane deformation with- 
out in-plane shear. The remaining deformations do not 
couple to fluid vorticity so we may set \ = hereafter. 
We expand the two remaining functions p, <E> in terms of 
solid spherical harmonics in the interior and exterior of 
the spherical membrane: 



p in = 


£*&(*) 

im 


Y tm {9,<f>) 


(29) 


$ in = 


im 


Y im (e,4>) 


(30) 


pout _ 


im 




(31) 


<jjOUt _ 


£*A?(*K 


-'- x y< m (M). 


(32) 



In terms of the basis of solid spherical harmonics we may 
write the velocity field in the exterior and interior of the 
sphere as a sum over angular modes n = (£, m) 



v = £ t V x ( r x«) + 
(n + 3) 



(33) 



n— — oc 



2rj(n + l)(2n + 3) 



?7(n + l)(2n + 3) 



Using Eqs. and |[2"9" ]) -([3"2 J) we find that the continuity 
of the radial fluid velocity implies 



£R l+1 
2 V (2£ + 3) 



pin 



L ff^in , (^+1)^ Pirn 

R hn + 2r/(l - 21) RZ+ 1 



R 



0. (34) 



The tangential velocity of the fluid in the absence of vor- 
ticity, on the other hand, is given by 



(35) 



Cm 



2 V (1 + 1)(2£ + 3) 



and 



,out 



(36) 



(2 - £)r 2 
2rj£{\ - 2£) 



,-i-i 



where we have used Eqs. ([29 j) -(|33 jl . 

We now expand the membrane deformations in spher- 
ical harmonics in a manner analogous to our expansion 
of the fluid velocity field by writing 



u>(M,t) = £u>/m(*)5WM) (37) 

im 

0, t) = £ ** ro (f)iW0, 0). (38) 



Using Eqs. ([37]) and fl38j) along with Eqs. ((351) and |[36J) 
in Eq. (|23|) . we find that the no-slip condition on the tan- 
gential velocity field at the membrane gives two equations 



and 



(I + 3)R e + 2 
2 v (£+l)(2£ + 3) 



P 1 







(2 - £)i? 2 



pout 

r im 



R^ 1 2f]£(l - 21) RC+ 1 



Im 



0. 



(39) 



(40) 



where the ' denotes the time derivative. 

We have now accounted for the kinematic boundary 
conditions on the fluid and membrane velocity at the sur- 
face of the sphere. We now require the balance of radial 
hydrodynamic and viscoelastic stresses at this surface to 
satisfy Darcy's law. To ensure force balance at the mem- 
brane we demand the equality of the vectors in Eq. (|25[) . 
Clearly this results in three scalar equations, but, it is 
computationally convenient to chose those scalars in a 
rather nontransparent manner [25| . We first set the com- 
ponents of the forces normal to the membrane equal to 
each other: 



f F 



fluid 



= f • F. 



(41) 



We then set the divergence of the two force fields equal 
to each other at the surface of the membrane using the 
vector identity 



(r-V) 



'fluid 



rV-F 



fluid 



-RV-F. 



r=R 



(42) 

Finally we also set the curls of these vector fields equal 
to each other on the surface of the membrane, 



[r ■ V x F fluid ] r=fl = r ■ V 



(43) 



These curls, in fact, vanish since we have excluded in- 
plane shear deformations of the membrane and vorticity 
of the fluid so this additional boundary condition is au- 
tomatically satisfied. 

We now enforce Eq. (14"T]) . The radial component of the 
forces on the membrane due to the fluid [25j | is given by 



r • Pfluid = V 



Cm 



-2{i +!)(£ + 2) $gjf(t) 



R 2 RZ+ 1 
31-1 P%£(t) , 2^-1) 



ry(2£-l) 
£ 2 - £ - 3 



RC+ 1 



R 2 



R^Tm(t) 



?7(2£ + 3) 



R e Pt(t) 



Cm { 



(44) 



The normal component of the viscoelastic forces on the 
membrane can be written in terms of w and VP, the radial 
and in-plane longitudinal deformations of the membrane. 
These forces are computed by a derivative of T\ : 



F r 
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r = r 
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- W (2/R-co)- W v, 



2K 



- k(A ± + 2/R 2 ) w-— A ± * 
R 



(45) 
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For the case of a viscoelastic membrane Eq. (|45j) involves 
an integral over the deformation history of the mem- 
brane. In the frequency domain, however, we incorporate 
the viscoelastic case by letting the elastic moduli become 
frequency-dependent complex quantities. The same ap- 
plies to the Eq. (|49|) below. The spontaneous curvature 
Co appearing in the force law Eq. (|45|) generates an ad- 
ditional, deformation-independent normal force on the 
membrane that is compensated by internal stresses in 
the membrane. Since we do not consider the nonlinear re- 
sponse of the elastic material such internal stresses do not 
affect the dynamics. Consequently, we set 2/R — cq = 
in the remainder of this work. 

Expanding Eq. (j45f in spherical harmonics and equat- 
ing these normal viscoelastic forces to the hydrodynamic 
ones from Eq. (|44|) as required by Eq. (|4ip leads to 



+ (£p+»'0-u' + 4) 



f- 

\R 4 

e+u-i 

(1 - 21) #+i 



R 2 

2 v (£+l)(£ 



2) n% 



{21 + 3) 



2*7 



R 2 
R 2 



R^Tm = 0. (46) 



Using the same analysis to demand force balance for 
the tangential components of the forces, we determine 
LHS of Eq. gU) to be 



[(r-V)(--F flnid ) - rV-F fl 
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RW 



r](2£ + 3) 



(47) 



The RHS of Eq. (|47|) can be written in terms of the de- 
formation field of the membrane, 



f x (f x F) — 



Vi 



2K 

IT 



SJ 7 ! _ 



2/i 
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(48) 



and 



-R V • F = ^-r + 2k [V 2 + 4t ) w- 2KV 2 , w 



R 2 



R 2 



fl(i^ + /i)Vi*-%Vi* + ^-Vi*. (49) 
R 



~R 



Expanding Eq. (|49[) in spherical harmonics and equating 
it to Eq. E3I we find 



^ 3 + 2^ + 6 , in 2r?l(l-l) 2 



(21 + 3) 



R 2 



£ 3 + 3l 2 + 51 - 3 P?£ 2t?(1 + 1)(1 + 2) 2 ggg 
(1 - 2£) + R 2 Rt+i 

(2 A {£ + 2f(£ - I) 2 + 2 J + 1) + 8 J) «* m 

+ -^(V + m)^+1) 2 -2m^+1) 

+ 4i4T£(£+l))^ ro = 0. (50) 

We use Darcy's law to set the normal velocity dif- 
ference between the fluid and the membrane, Eq. (|27[) . 
This generates an ordinary first-order differential equa- 
tion. We look for solutions having a time dependence of 
the form ~ e~ lu)tmt . Using this assumption in Eq. f2?]). 
we find our sixth and the last relation between the dy- 
namical variables, 
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(51) 



We now have a set of six algebraic and differen- 
tial equations governing the six dynamical variables 
(#, w, $ in , P in , $ out , P out ). These six variables are gov- 
erned by Eqs. (JSJ), ((39]), gDJ), I0BJ), © and (51). In this 
a damped system, we expect to find solutions having a 
negative imaginary part of W£ m , which we will interpret 
as the decay rate of a mode describing the coupled dy- 
namics of the spherical membrane and fluid system in 
the case of vanishing in-plane membrane shear and fluid 
vorticity. 

Since we want to focus on the dynamics of the mem- 
brane, we use Eqs. ([34"]). ([39| . ([40]) . |46|) to eliminate the 
four variables p, $ associated with the fluid dynamics in- 
side and outside the membrane. As a result we find a 
system of two algebraic equations for wg m and $f m de- 
scribing the compression and bending of the membrane 
respectively: 



.4 



Win 



= 0. 



(52) 



The a>£ m -dependent components of the A matrix are 
listed in Appendix [C] The condition for the existence 
of nontrivial solutions to Eqs. (f52"| determines the decay 
rates of the deformation modes of the membrane coupled 
to the surrounding fluid. 

To discuss these decay rates of the membrane defor- 
mation modes, it is convenient to introduce the following 
dimensionless variables 



UJ£T]R 3 _ II R 2 
, A 1 = : 



aR 
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KR 2 



(53) 
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D. The modes of the shell 

In Fig. [2] we plot the decay time Tp — \ jiur}' of the 
first mode as a function of the order of the spherical har- 
monic £ > 1. To demonstrate the effect of permeation on 
these decay rates we show three curves corresponding to a 
highly permeable membrane (a = 0.1), one of intermedi- 
ate permeability (a = 25) and a completely impermeable 
one (a — oo). All three cases correspond to a perfectly 




FIG. 2: (color online) Dimensionless relaxation times of the 
Lennon Brochard mode fi for an elastic membrane vs the 
mode number £ for various permeation parameters a. The 
dimensionless elastic parameters of the membrane are p, = 75 
and K = 225, which corresponds to the Poisson ratio a = 0.5 
and the radius to thickness ratio R/h = 5. =0 for a — 
for all I 

elastic membrane so that p, and K are real and frequency 
independent. For £ ^S> 1 the characteristic wavelength of 
the deformation is much smaller than the radius of cur- 
vature of the sphere so we should expect to recover the 
flat membrane result of Lennon and Brochard [26| . In- 
deed, for large £, we find that, t ( > ~ £ so that the 

decay rate (f^l ~ ^ ~ (x) 3 ~ 9 s wnere ^ * s tne 
wavelength of the mode. As expected based on our anal- 
ysis of the plane permeable membrane, the transition to 
the Lennon and Brochard behavior is more apparent for 
the impermeable case (a = oo) and becomes less distinct 
as a decreases. It is not possible to consider the case 
a = when discussing the relaxational dynamics of the 
LB mode as the decay rate of this mode becomes infinite 
as a — > for all £. 

We parameterize the relative degree of bending versus 
compression associated with these modes by introducing 
their ratio 

w tm /9 tm = Z {p )(£), (54) 

where p = 1,2 indexes the two linearly independent 
modes of the system. We plot in Fig. [3]this ratio Z^\£) 
for the first mode whose corresponding decay times are 
shown in Fig. [2] . Examining this plot it is clear that this 
mode is, indeed, the analogue of the Lennon-Brochard 



bending mode since \Z^ l \£)\ 3> 1. This identification 
of the first mode with bending becomes more exact at 
higher £ where the curvature of the membrane becomes 
less significant for the dynamics. The overall negative 
sign for this result is due to the fact that where the mem- 
brane bends outward the material expands and where it 
curves inward the material compresses. Due to the over- 
whelming bending character of this mode we will refer to 
it as the Lennon-Brochard (LB) mode. 




50 100 150 200 

t 

FIG. 3: (color online) Ratio of bending to stretching Z^(€) 
in the Lennon-Brochard mode vs. mode number I of a per- 
fectly elastic membrane. The membrane elasticity parameters 
are identical to those used in Fig. [2] The values of the dimen- 
sionless permeation coefficient are listed in the figure legend. 

We now turn to the second mode and plot in Fig. 2] its 
dimensionless decay time f^ 2 ' = 1/iuj^ as a function of 
£. We use the same elastic parameters p, and K, as in the 
plot of the decay times of the LB mode and show curves 
corresponding to an impermeable case (a = oo), one of 
intermediate permeability (a = 25), and a case where 
the fluid passes through the membrane without any re- 
sistance (a — 0). When £ ^> 1 all these cases converge 
onto a single curve so that rf^ ~ £~ 1 independent of a. 

As fg ~ £~ 1 for large £, this mode cannot be con- 
trolled by membrane bending at high wavenumber. Turn- 
ing to Z^ 2 \£) - see Eq. (|54)l -, which we plot in Fig.[5l we 
see that this mode is dominated by the in-plane compres- 
sion of the membrane. In fact, this mode is the remnant 
of longitudinal sound (LS) in the elastic membrane that 
is over-damped due to the coupling to the surrounding 
solvent. It is interesting to note that only at large £ the 
(LB) mode becomes strongly bending dominated while 
the longitudinal sound mode (LS) only loses its bending 
character at large £. It is well known that bending and 
longitudinal sound are orthogonal normal modes of the 
linearized membrane/solvent system in the case of flat 
membrane 27]. Examining the results for the sphere we 
find that geometry couples these deformations, but, at 
high wave number (or large £) where the wavelength of 
the deformation is much less than the radius of curvature 
R, the two orthogonal modes of the sphere approach the 
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FIG. 4: (color online) Dimensionless relaxation time ff^ in 
the longitudinal sound mode vs. the mode number I for purely 
elastic membranes having various dimensionless permeation 
parameters a. The dimensionless elastic parameters of the 
membrane are \x = 75 and K = 225 as in Fig. [2] 




100 



FIG. 5: (color online) Ratio of bending to stretching (£) in 
the longitudinal sound mode vs. mode number I of a perfectly 
elastic membrane. The membrane elasticity parameters are 
identical to used in Fig. [3] The values of the dimensionless 
permeation coefficient are listed in the figure legend. 



k + R 2 K 
R 1 



(56) 



The dimensionless relaxation rate of the t = LB mode 
is then given by Eqs. (|55|) and (|56|) . We find 



f (2) 



1 



a 



4 1 + K 



(57) 



It is clear that the decay rate of this mode vanishes in 
the limit of an impermeable (a = oo) sphere as volume 
changing deformations are not allowed. As the sphere be- 
comes more permeable the decay rate increases becoming 
infinite in the limit of a "ghost membrane" (a — 0) since 
we have neglected the membrane's inertia. 

To compare these results with experiments it is useful 
to recognize that our two-dimensional solid membrane 
is an idealization of a thin sheet of thickness h whose 
material properties can be defined in terms of three- 
dimensional Lame coefficients /i3, A3. Using standard 
results [U HtJ we may then reexpress our dimensionless 
elastic parameters in terms of the geometry of the system 
and the Poisson ratio of the membrane material a. We 
find 



/1 = 6 



R 



R 



(I -a) 



while the bending modulus is given by 
/x 3 /i 3 A 3 + ^3 



3 A3 



2^3 



(58) 
(59) 

(60) 



We return to this point with regard to the microrheology 
of viscoelastic spherical membranes in section IIIII 



III. THE NANOINDENTATION RESPONSE 
FUNCTION 



simple mode structure of the flat membrane. 

Before closing our discussion of the mode structure of 
the spherical elastic membrane, we must consider the 
I = mode of the membrane, which corresponds to 
purely radial motion of its surface. This mode must 
be handled independently of the others because of the 
incompressibility of the solvent. Because of that incom- 
pressibility, there can be no corresponding radial motion 
of the solvent so that the membrane relaxes against a 
background of quiescent fluid. In this mode the radial 
deformation obeys 



Using the above analysis of the mechanics of porous, 
viscoelastic spherical shells, we now determine the mi- 
crorheological signature of viscoelastic, porous shells. As 
discussed above, these results are relevant to the study of 
the mechanics of viral capsids, colloidosomes, and porous 
lipid bilayers such as those that contain membrane-bound 
pore-forming proteins. 

In the active AFM-based mechanical measurement, the 
shell is deformed by antipodally applied forces directed 
radially toward the center of the shell. To measure the 
mechanics of the shell the distance between the point 
of force application and its antipode is measured. This 
distance, due to the application of a sinusoidally varying 
force Foe~ iuJt , will deviate from the shell diameter with 
an amplitude 



with no angular dependence, and with an effective spring 



D(uj) = 2R-H{w)F , 



(61) 
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which defines the response function H(ui). In the passive 
microrheological experiment, on the other hand, the ther- 
mally excited fluctuations of 2R — D give the same infor- 
mation via the fluctuation-dissipation theorem (38l . |39| . 

To calculate H(uS) we apply equal and opposite forces 
to the "north" and "south" poles of the sphere, i.e. at 
polar angles 9 = 0, ir. The force balance boundary con- 
dition Eq. (14"T]) then becomes 

r • Ffluid + wi (w) + 7T 2 (w) = f • F, (62) 
where are the external normal stresses on the surface 

7Ti(o;) = -<J(0) e - iwt rF„/ fl 2 (63) 
tt 2 (w) = < 5(6'-7r)e-^ t fF /a 2 . (64) 

As shown above we will consider point forces so that 
the area of the sphere contacted by the AFM tip a — > 
at constant Fq/o 2 — ttq. Using the completeness of 
spherical harmonics on the unit sphere we express these 
externally applied stresses in terms of spherical harmonic 
modes: 

7r(t, 9, 0) = -r tt Ko(0, 0) + F/ (tt, 0)] J$ m (0, fle"*"* 

(65) 

Thus the AFM-induced pinching of the shell couples to 
all the spherical harmonic modes of the system discussed 
above. The effect of the finite-size of the AFM tip can 
be accounted for by putting in a large--^ cutoff in the 
sum shown in Eq. Q65[) . For reasonably small tip, this 
effect must be negligible. Using these previous results 
we compute the radial deformation amplitude associated 
with each I— mode at the north pole of the sphere. In 
terms of the dimensionless frequency, Q = uorjR z / k, this 
radial displacement of the sphere at the north pole is 
given by 

w t (Q) = 2— tto[1 + (-1)"]17 (0, 0)6(0), (66) 

where Q{Co) is the frequency-dependent complex radial 
compliance of the shell due to a radially applied normal 
stress distributed over the sphere as Yi (9). This com- 
pliance results from both the Lennon-Brochard (LB)and 
longitudinal sound (LS) modes of the shell discussed 
above. The complex, frequency-dependent response 
function can be expressed in terms of five dimensionless 
functions of the angular momentum of the mode £: 



k x {l) + iuk 2 (£) 



(67) 



where due to their length, these functions Ai(£), A$(£) 
are reported in Appendix [Dl As shown there, it is 
through these functions that the ^—dependent response 
function acquires its dependence on the dimensionless 
permeation parameter a and on the elastic constants of 
the membrane written in terms of the ratio R/h and the 
Poisson ratio a as shown in Eqs. I[58p and (|59[) . 



Summing over the I— modes using Eqs. ((65)) and (|66|) 
we find the radial displacement of the shell at (0, 4>) to 
be 

W (u>, 9^) = ^^J2 Y ^ 9 ' + (-±YM")y; (o, 0). 

K * — ' 

e=o 

(68) 

Thus the amplitude of oscillations of the diameter is given 
by 

D(Q) = 2R + w(u>, 0,0) + w(Q,n,0). (69) 

Using the definition in Eq. flST} and Eqs. ([57J>. ([55)). and 
(1691 ), we solve for the response function in terms of the 
point force response of the shell. We find that the diam- 
eter response function is given by 

2R 4 00 

H(w) = £[1 + (-l)%(d,)\Y m (0, 0)| 2 . (70) 

7=t 



A. An elastic shell 

We plot in Fig. [HUa) the real (in phase) and imaginary 
(out of phase) parts of the diameter response function 
H(Cu) for purely clastic shells of varying porosity. In an- 
ticipation of future fluctuation-based microrheology ex- 
periments we plot in Fig. 5(b) the predicted power spectra 
for the thermal fluctuations of this variable, computed 
using the fluctuation-dissipation theorem [12j. Examin- 
ing the figure we note that the response function can be 
characterized as having a low-frequency elastic plateau 
followed by a cross-over to a viscous-dominated decay 
at high frequencies. As the dimensionless permeation 
coefficient is varied from that of an impermeable shell 
(a = oo ) to one that allows fluid permeation with no re- 
sistance (q = 0), the dimensionless cross-over frequency 
moves from rs 10 to ss 10 2 . We return to predicted val- 
ues of this cross-over frequency for various physical sys- 
tems in section IIVI Within the low-frequency plateau 
where much, if not all, the dynamical data is likely to 
be taken, the main effect of permeation is to shift the 
value of the response function by approximately a factor 
of two. Thus, in order to correctly extract the bending 
and compression moduli of the shell to better than this 
factor of two, an accurate account of the permeation must 
be made. In the passive microrheology experiment, the 
effect of permeation on the thermal power spectrum of 
diameter fluctuations is more pronounced. Ignoring per- 
meation in this sort of measurement leads to quantitative 
inaccuracies on the order of 4 as shown in Fig. \&[b) . We 
now consider the case of a viscous shell. 



B. A viscous shell 

The case of a viscous shell can be physically realized 
with giant unilamellar vesicles (GUVs) that contain pore- 
forming transmembrane proteins or protein complexes. 
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FIG. 6: (color online) The response function of a porous, elas- 
tic shell. In panel (a) we plot the real (Re[H(uu)], solid lines) 
and imaginary (Im[H(u))], dashed lines) parts of the response 
function vs dimensionless frequency Q for various values of the 
permeation parameter a. The remaining material parameters 
are: R/h = 5, Poisson ratio a — 0.5. In panel (b) we show 
the thermal power spectrum (Im[H(p)]/{j) using the same 
material parameters. 



This system allows for the largest variation of permeation 
coefficients of any that we study. By varying the num- 
ber of such proteins in these protein/lipid complexes, one 
may vary the porosity parameter from essentially a — oo 
(no pore-forming protein) to finite values. Using the solu- 
tion for the permeation coefficient found in Appendix [XI 
it should be possible to reach values of a = 10 for the 
lipid vesicle systems. 

To determine the dynamics of a purely viscous shell we 
replace real, frequency- independent shear modulus of the 
shell by a viscous response: [i — » — iw7y m , with rj m being 
the two-dimensional surface viscosity of the the mem- 
brane or bilayer. To allow any radial deformations of the 
impermeable bilayer, there must be some area change 
of the surface. One may imagine that the surface of 
a tensed GUV is essentially inextensible, however, even 
such tensed structures contain a reservoir of extra area 
hidden in the (small-scale) thermal undulations of the 
surface [13 . [40| . Thus the two-dimensional bulk mod- 
ulus of the viscous shell is dominated by a large, but 
finite elastic component. We expect that there may well 
be a dissipative or viscous response of the shell to area- 
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FIG. 7: (color online) Response function of a GUV contain- 
ing pore-forming protein complexes. In panels (a) and (b) we 
show the real (Re[H(ui)]) and imaginary (Im[i?(a))]) parts of 
the response function vs. dimensionless frequency uj respec- 
tively using various values of the permeation parameter a. 
The other dimensionless parameters were:R/h — 20, a — 0.5 
rim/Rri — 0.1/io, jJ-0 and K defined in Eqs. (|58[l and (159[) . 
For comparison, black solid line shows the real part of the 
response function for an elastic shell with R/h — 20, a — 0.5, 
a — 50. The black dashed line corresponds to the imagi- 
nary part of the response function of an elastic shell with the 
same material parameters as in panel (a). In (c) we plot the 
thermal power spectrum {lm[H{uj)\/uj) as a function of the di- 
mensionless frequency u> using the same material parameters 
as in panel (a). Black line - elastic shell. 



changing deformations, but we ignore these subdominant 
corrections. In Fig. [7{a), we plot the real and imaginary 
parts of the response function for a porous (a = 50), 
purely viscous shell with a finite area compression mod- 
ulus, K. In part (b) of the same figure, we calculate the 
power spectrum of thermal fluctuations of the diameter 
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of such an object. 

The low frequency response of the porous, viscous 
membrane exhibits an elastic plateau due to its elas- 
tic response to both bending and compression. The 
viscous shell, however, admits a new intermediate fre- 
quency response regime associated with viscous dissipa- 
tion within the fluid membrane. In this intermediate 
frequency regime, the real part of the response function 
decays as d>~ - 66 over about two decades in frequency. 
Examining these results we note that there are dramatic 
qualitative differences between the response functions of 
viscous and elastic shells at least within this intermediate 
frequency range. 

In Fig.[7]Jb) we see the low- frequency peak in the imag- 
inary part of the response function associated with these 
internal (to the membrane) shear modes. At still higher 
frequencies, both the real and imaginary parts of the re- 
sponse function are dominated by dissipative stresses due 
to solvent flow and permeation through the pores. This 
can be most easily seen in the correspondence of the high- 
frequency peak of the imaginary part of response func- 
tion of the viscous membrane with only peak in the cor- 
responding part of the response function for the purely 
elastic shell. In this high frequency limit a viscous mem- 
brane and an elastic shell show essentially an identical 
mechanical response. This is reasonable since the high- 
frequency dynamics of both systems are dominated by 
solvent flow and permeation. 

Passive microrheological measurements at low frequen- 
cies can easily distinguish between a viscous and an elas- 
tic shell as can be see by the difference between their 
power spectra shown in Fig. 0c). The effect of perme- 
ation, as can be seen by comparing the curves, is prin- 
cipally in shifting the hight of the plateau region and 
slightly altering the transition frequency to the solvent- 
dominated, high-frequency regime. 



C. Viscoelastic shells 

To explore the role of viscoelasticity in the mechanics 
of a porous spherical shell, we consider two archetypal 
cases. We first calculate the response function for a vis- 
coelastic shell where the shear modulus has a single re- 
laxation time, t s . This Maxwell model [23[ of the shear 
modulus can be written as 



MM = Mo 



1 — iU!T s ' 



(71) 



where f s is the dimensionless shear stress relaxation time. 
Deformations occurring at a frequency scale much lower 
than the inverse stress relaxation time relax viscously. 
The high frequency dynamics, on the other hand, see an 
elastic shear response of the material with shear modulus 
Mo- 

In Figs. [Eta) and (b) we plot the real and imaginary 
parts of the response function respectively for a viscoelas- 
tic shell having a complex, frequency-dependent shear 



modulus of the Maxwell form. We take the shear stress 
relaxation time f s to be within the range of decay times 
of the various £— modes of the system. Thus, the small 
£— modes of the shell have long decay times relative to 
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FIG. 8: (color online) The response function of a shell with 
a viscoelastic shear modulus. In panels (a) and (b) we show 
the real (Ke[H(u)]) and imaginary parts of the response func- 
tion vs. dimensionless frequency uo respectively using various 
values of the permeation parameters a. The other dimen- 
sionless parameters were: R/h — 20, a = 0.5, f 3 = 0.1, p,o 
and K defined in Eqs. (|58[) and ((59}. For comparison, the 
black lines shows the real (solid) and imaginary (dashed) parts 
of the response function for an elastic shell with R/h — 20, 
a = 0.5, a — 50. In (c) we plot the thermal power spectrum 
(Im[H(u))]/ui) as a function of the dimensionless frequency Q 
using the same material parameters as in panel (a). Black 
line - purely elastic shell. 

the stress relaxation time and thus experience viscous 
relaxation dynamics, where the higher £— modes of the 
system decay fast enough to be affected by the elastic 
shear response of the system. 

It is that for a viscoelastic shell the intermediate fre- 
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quency response regime associated with the viscous sys- 
tem remains, but the frequency dependence of the real 
part of the response function is weakened: ~ uj~ 0A4 . This 
appears reasonable as the viscoelastic case effectively in- 
terpolates between the response of an elastic shell (high 
^-modes) and a viscous one (low ^-modes). The slope 
of the real part of the response function (see Fig. [51(a)) 
in the intermediate response regime now is a-dependent. 
Thus, to use the response function to extract properties 
of the viscoelastic response of the shell is difficult without 
an accurate model of solvent permeation through it. 

In Fig. [5(c) we plot the thermal power spectrum for 
this model. As long as the shear stress relaxation rate is 
small compared to the cross-over to the high-frequency 
terminal behavior of the system, it is indeed possible to 
observe the viscoelastic response of the shell by exam- 
ining the thermal power spectrum. This identification, 
however, relies on a separation of time scales so that 
f s remains long compared to the cross-over timescale to 
viscous-dominated behavior. As we discuss in section llVl 
this separation of time scales should be easy to achieve 
in many physical systems. 

Finally, we study perhaps the most potentially inter- 
esting example of a viscoelastic shell on the nanoscale. 
In some viral capsids the individual capsomeres can un- 
dergo a conformational change that changes their effec- 
tive cross-sectional area. These internal degrees of free- 
dom related to volume-changing allosteric transitions of 
the proteins making up the shell of the virus allow for new 
dissipation associated with area-changing deformations. 
To model this system we include a complex frequency- 
dependent two-dimensional bulk modulus of the form 



1 



-ILUTB 

I 

1 — iuJTB 



(72) 



where Kq sets the overall modulus scale, tb is the inverse 
rate of conformational change, and gKo sets the maxi- 
mum of the dissipative part of the modulus at uitb — 1- 
The typical time scale for allosteric transition in proteins 
is on the order of microseconds to milliseconds. The ra- 
tio of the zero-frequency elastic part of the area modulus 
Kq to the dissipative contribution at the frequency of 
these allosteric transitions is controlled by the constant 
g. Since we know of no quantitaive measurement of these 
parameters, we sweep the value of g from zero, which cor- 
responds to a purely elastic shell as shown in Fig. [SJ to 
g=l, where this dissipative contribution to the stress 
dominates the zero-frequency elastic contribution to the 
area modulus. It is clear from Fig. [9] that even a small 
dissipative response to compressive stresses at time scales 
associated with allosteric transitions leads to a clear and 
measurable change in the response function of the shell. 



IV. SUMMARY AND DISCUSSION 

We have developed a theoretical model for the dy- 
namics of a porous spherical shell immersed in a vis- 
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FIG. 9: (color online) A simple mechanical model of a vi- 
ral capsid in which the capsomers undergo a cross-sectional 
area changing allosteric transition. In panel (a) we plot the 
real and imaginary parts of the response function for different 
Maxwell solid models of the shell - see Eq. (|72|) . In panel (b) 
we plot the expected thermal power spectrum (lm[H(Q)]/uj). 
In all cases the material parameters of the shell are given by: 



R/h = 5,t b = 10 _;i s, Y = 140MPa, a = 0.5, 
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cous solvent and used it to examine the utility of 
nanoindentation-based rheological measurements. In the 
limit of a flat, porous membrane immersed in a vis- 
cous solvent, we calculated the relaxational dynamics of 
the bending or undulatory modes of the surface. These 
modes are well-understood in the limit of a impermeable 
surface. In that limiting case, the long-range hydrody- 
namic flows introduce a nonlocal coupling of the stress 
and displacement along the surface. It has long been 
known that the hydrodynamic coupling shifts the dis- 
persion relation of the undulation modes of wavevector 
q from their expected "local drag" form of uj ~ — iq 4 
to correct answer of uj ~ — iq 3 incorporating the hydro- 
dynamic cooperativity of the system as first proposed 
by Lennon and Brochard For the case of perme- 

able membrane, the dispersion relation of the undulatory 
modes smoothly crosses over from the local drag form to 
the usual Lennon-Brochard solution as the wave num- 
ber is varied. The cross-over length is set by £ ~ rj/a. 
At lengths much larger than i the solvent permeation 
generates a subdominant correction to the usual dynam- 
ics of the undulatory modes of a membrane. Below this 
length, however, permeation effectively destroys the hy- 
drodynamic coupling along the membrane leading to a 
new dynamics consistent with the naive, local-drag ap- 
proximation. 



14 



In order to consider the mechanics of spherical perme- 
able membranes (of radius R), we examined both com- 
pression and bending deformations of the material since 
curvature couples these linearly independent modes of 
the flat membrane. The hydrodynamic interactions of 
the membrane are also less simple in this geometry. Using 
a solution for the solvent dynamics both inside and out- 
side the spherical shell, we determined the over-damped 
normal modes of the system. Owing to the spherical 
symmetry of the problem, these modes are proportional 
to the spherical harmonics and thus can be indexed by 
the usual angular momentum variables, (£,m). At small 
angular momentum, i.e. I ~ 1, the two linearly inde- 
pendent modes of the membrane each combine an undu- 
latory and compressional character. At higher angular 
momentum corresponding to higher wavevectors on the 
surface q ~ l/R, the dynamics approach those of the flat 
interface where these two modes acquire a dominantly 
undulatory or compressional character respectively. Us- 
ing this short length scale behavior, we label the modes 
as either the Lennon-Brochard (i.e. undulatory) or com- 
pressional type. 

We use these dynamical calculations to study the ex- 
pected finite-frequency response of a variety of essen- 
tially spherical, viscous, elastic, and viscoelastic porous 
shells. Such systems include colloidosomes, nanoparti- 
cle networks, vesicles containing membrane-bound pore- 
forming protein complexes, and viruses. The principal 
questions that we address are how does porosity affect 
the mechanical response of viscous, elastic, and viscoelas- 
tic spherical shells and specifically, for the case of viral 
capsids, is it possible to mechanically detect allosteric 
transitions in the capsomeres making up the protcinous 
outer shell of a virus? 

To address these questions, we need only determine the 
characteristic frequencies associated with the nondimen- 
sionalized curves shown in Figs. El El El ancl 13 A few 
general comments are in order. First, the characteristic 
frequency associated with the transition between the low- 
frequency plateau and high-frequency, fluid dominated 
decay of the response function is typically large enough 
in viruses to render the high-frequency region unobserv- 
ablc using current methods. Expressing the frequencies 
/ in dimensional form using the dimensionless angular 
frequency Q and the relation 



Y 



24tt?7(1 - a 2 



{h/Rf, 



(73) 



we find that the characteristic cross-over frequency for 
a CCMV virus 11] to be on the order of 10 8 Hz. For a 
larger and mechanically softer liposome, however, this 
transition frequency is ~ 100kHz. For colloidosomes 
and other synthetic structures on the scale of one to 
ten microns we estimate this transition frequency to be 
~ 10 2 — 10 4 Hz. The general trends are simple to under- 
stand from Eq. (|73|) . Larger (large R) and thinner- walled 
(small h) spheres have a lower transition frequency, while 
those spheres made of less elastically compliant materials 



(i.e. higher Young's modulus Y) have higher transition 
frequencies. 

We expect, based on the above estimates, that the 
transition to a permeation-dominated response is not ob- 
servable in viral capsids. One advantage of this point is 
that an observed increase in the imaginary (out-of-phase) 
response function for viral capsids implies the existence 
of internal dissipative modes of the structure that are 
likely to be related to allosteric transitions of capsomeres 
or in their rearrangement relative to one another. We 
show in Figs. 18191 two possible rheological signatures of 
such internal dissipative modes. In the former we sup- 
pose that there is a large elastic response of the shell to 
area-changing deformations, but that the shear response 
of the material is viscoelastic. One might expect such 
a mechanical response if the individual capsomeres irre- 
versibly rearrange in response to applied stress. In the 
latter case we explore a system having a Maxwell solid 
response to area compression, but an elastic response to 
shear. Such a material response might be expected if the 
individual capsomeres undergo allosteric transitions that 
change their effective cross-sectional area. In Fig. HI we 
see that the Maxwell shear response leads to measurable 
corrections to the rheological spectrum of the object. 

The AFM-based mechanical probe is also highly sen- 
sitive to a viscoelastic response of the area compression 
modulus as might be due to allosteric transitions in the 
capsomeres. The expected separation of time scales be- 
tween the lower frequency allosteric transitions (10 3 Hz 
to 10 6 Hz) and transition to permeation dominated dy- 
namics (10 8 Hz) demonstrates that this measurement is 
well-suited to exploring such capsomere conformational 
changes. 

The permeation parameter is not irrelevant for making 
quantitative measurements. This parameter shifts the 
value of the elastically dominated low-frequency plateau 
of elastic shells (such as viruses) by as much as a fac- 
tor of four. This dynamical response function is princi- 
pally controlled by the bending modulus of the shell so 
that AFM-based measurements of that quantity can be 
in error by as much as a factor of ten if permeation is 
neglected in the interpretation of the data. The effect of 
solvent permeation is generically more important in sys- 
tems that are both larger and mechanically more com- 
pliant. Thus the effect of solvent permeation can play 
a large role in the mechanical response of large GUVs 
containing pore-forming proteins. 

The combined effects of solvent permeation and vis- 
coelasticity remains to be fully explored in a wide variety 
of biological and synthetic nanoscale shells and vesicles. 
Based on our present work, accounting for solvent perme- 
ation will be important for the quantitative interpreta- 
tion of such rheological data. By taking such effects into 
account it may be possible to make sensitive measure- 
ments of the collective effects of protein conformational 
change in viral capsids and to probe the mechanics of a 
number of other similar structures. 
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APPENDIX A: THE PERMEATION 
COEFFICIENT 

To determine the permeation coefficient we treat the 
porous membrane of thickness t as having an array of 
pores or radius rg- We use the well-known result for the 
fluid flux Q through such a tube of length t and radius 
ro due to a pressure difference AP across the tube to 
write (3 



nrf AP 
8r/ t 



(Al) 



We then suppose that the number of such tubes is k per 
unit area of the membrane and assume that the total flux 
through these tubes is simply additive. Then the total 
flux through the area L 2 is thus given by 



total flux = kL 2 Q, 



(A2) 



and the average velocity of fluid coming out of the mem- 
brane is 



fc- 



■nri AP 



8?7 t 

On the other hand, we assert from Darcy's law that 

AP 



(A3) 



(A4) 



From Eqs. (|A3|) and (|A4|) we find the permeation coeffi- 
cient to be given by 



a 



8rjt 



(A5) 



Expected values of the permeation coefficient for a va- 
riety of micron and nanoscale hollow, porous shells are 
given in Table 1. In the text we introduce a dimension- 
less form of the permeation coefficient a that is defined 
by 



a 



a ■ R 
V 



\Rt 



where R is radius of the sphere and r\ the viscosity of the 
surrounding solvent. 



APPENDIX B: UNDULATIONS OF A PLANE 
PERMEABLE MEMBRANE 

We calculate the fluid flow associated with the undu- 
lations of a permeable membrane. To begin we look for 



V-v = 

<9v 2 
p— = 77V v - Vp. 
at 



(Bl) 



driven by sinusoidal membrane undulations of the form 



iqx — iu)t 



(B2) 



These undulations apply an external normal stress on the 
fluid at the membrane located in the xy plane (z = 0) of 
the form 



iqx — iujt 



(B3) 



Correspondingly, we expect a solution for the fluid veloc- 
ity field v and pressure p of the form 



v = v (z)e i<?:E_i " t 
p = p{z)e iqx ~ lujt , 



(B4) 
(B5) 



but the fluid velocity and pressure fields must vanish far 
from the membrane z — > ±00 and satisfy at the mem- 
brane 

»+(0)=«-(0)=0 (B6) 
v+(0)=v-(0) (B7) 

Additionally, the fluid incompressibility condition gives 

dv z 



dz 



+ v x (z)iq = 0. 



Using these in Eq. (|B2[) equation we find 



-piuv x (z) =r)[-q v x + 



2 d 2 v. 1 



dz 2 



(B8) 



iqp(z) (B9) 



dz 



The general form of the solutions to these ordinary dif- 
ferential equations is 



v x (z) = ae Xz v z {z) = be Xz p{z) = ce* 



(Bll) 



where the incompressibility condition Eqs. (|B1[) and (|B8[) 
requires that the undetermined constants are related by 



(A6) w hile Eqs (|B9]) and (|B10|) give four solutions 



\2 2 2 

A =q ,q - 1 — . 

V 



(B12) 



(B13) 



Using the boundary conditions we find solutions of the 
form 



v ±(z) = T iqA(e^ z - e Tlzyvx-i"t (B14) 
vf(z) = qA(e Tqz - l e Tizyq*-^t 
p ± (z) =TiupAe Tqz e iqx - iut , 



(B16) 
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with I = \J q 1 — ipuj/rj. 

Knowing the pressure and velocity fields, we require 
force balance at the membrane: 



7T = KV W, 



(B17) 



using the fluid stress tensors and membrane bending 
modulus discussed in the text. 

The permeability of the membrane sets the velocity 
difference between the fluid and the membrane at z = 0. 
Thus we also have 



a v z (0) - 



dw 

IT 



= (p -p 



(B18) 



These conditions set the two remaining unknown con- 
stants above, A from Eqs. (|B14|) - (|B16p and w q by requir- 
ing nontrivial solutions to the set of algebraic equations 



.4 



aq(l- - 



2 iujp 



+ w„ aitu = 



A 2iujp — Kq w q 



(B19) 



and 



A 21 =-3(2*+l)S 

+ 2(£-l) 2 (£ + 2) 2 (2£+l)- 



1 



tUJl 



-3(e-ifeu + 2) 2 -=- 

IUJ£ 



+ 12(2*+1) 



K 

iuj e 



+ 2£(£ + 2)(2£-\) 



aK 

iuji 



A 22 = i{2£ + l)(A£ 2 + 6£ - l)a 



U{l+l)(2£+V 



) K 



£ 2 {£+\){£ + 2){2£~l) — 

ILl>1 

(£- l)£(£ + 2)(2£ 2 + 3£ + 4) 



ap 



(C3) 



(C4) 



Working in the limit where fluid inertia is irrelevant ujp 
rjq 2 and solving for w q we find 



APPENDIX D: THE RESPONSE FUNCTION 



as discussed in the text. 



7T 


a 
2q V ' 


f 2 




2itoa — 


Kg 4 


a 
2qr 


+ 2 



Here we list the £ dependence of each function enter- 
(B20) ing into the diameter response function shown in Eq. (|67[) . 

These functions depend also on the dimensionless perme- 
ation parameter a, dimensionless shear modulus p, and 
dimensionless area compression modulus K. 



APPENDIX C: THE A MATRIX 

The eigenvectors and eigenvalues of the A matrix de- 
fined by Eq. (|52p determine the relaxation times and 
mode structure of the coupled membrane deformation 
and solvent flow modes. The components of this matrix 
are listed below: 



An = (2£ + l)(4£ 3 + 6£ 2 ~£ + 3)a 



(CI) 



(£ - lf(£ + 2) 2 {2£ + l)U£ 3 + 6t-£ + 3) — 

tut 



(£ - 1) 2 £(£ +!)(£ + 2) 2 (2£ + 1)- 



tUJ£ 



A(2£ + 1)U£ 3 + U 2 - £ + 3) (2£ + 1) — 

tut 

2£(£+l)(£ + 2) — , 
iu>i 



A 12 =£{£ + !) -3{2£+l)a 

+ 2 (2£+ 1){4:£ 3 + U 2 - £ + 3)- 



(C2) 



K 



3(£-i)l(£ + 2)^ 



Ki{£) = 2(1 + 2£){<o + 7£ + 7£ 2 )K 

+ £(£+ l){£ 2 + £ + A)aK 

+ (£-l)£(l + £){2+£)afL 
A 2 {£) = 2(1 + 2£)(3 + I + £ 2 )a 
A 3 (£) = (£-l) 2 £(l+£)(2 + £) 2 x 

(l + 2£)(3 + 2£ + 2£ 2 )K 

+ (£- l) 2 £ 2 (l + £) 2 (2 + £) 2 aK 

+ {£- 1) 3 (2 + £) 3 {l + 2f)(3 + 2£ + 2£ 2 )p 

+ (£-l) 3 £(l+£)(2 + £) 3 ap 

+ 4(1 - 1)(2 + £)(1 + 2f)(3 + 2£ + 2£ 2 )Kp 

+ A{£ - \)£{l + £){2 + £)aKp 

A 4 (£) = -(2£+ \)\{£- lf{£ + 2) 2 x 



(Dl) 



(D2) 
(D3) 



(D4) 



(21- l){2£+ l)(2£ + 3) 
+ 2(£-l) 2 £(£ + l){£ + 2) 2 a 
+ 4(2£-l)(2£+l){2£ + 3)K 
+ £{l + £){2£ 2 + 2£ - I) aK 

+ (£- 1)(2 + f)(3 + 2£ + 2£ 2 )ap 

A 5 (£) = ~(2£ - 1)(2£ + l) 2 (2£ + 3)a 



(D5) 
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